Small-scale turbulent dynamo 
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Kinematic dynamo theory is presented here for turbulent conductive fluids. We describe how 
inhomogeneous magnetic fluctuations are generated below the viscous scale of turbulence where the 
spatial smoothness of the velocity permits a systematic analysis of the Lagrangian path dynamics. 
We find analytically the moments and multi-point correlation functions of the magnetic field at 
small yet finite magnetic diffusivity. We show that the field is concentrated in long narrow strips 
and describe anomalous scalings and angular singularities of the multi-point correlation functions 
which are manifestations of field's intermittency. The growth rate of the magnetic field in a typical 
realization is found to be half the difference of two Lyapunov exponents of the same sign. 



It is believed that the magnetic fields of planets, stars 
and galaxies have their origin in dynamo action driven 
by motions of conducting fluids [|]-||. Inhomogeneous 
flow stretches magnetic lines amplifying the field while 
the field produces electric currents that dissipate energy 
and diffuse the field due to finite resistivity. The out- 
come of the competition between amplification and dif- 
fusion depends on the type of the fiow. We consider 
the long-standing problem of how turbulence excites in- 
homogeneous fluctuations of magnetic fleld |Q-^. Since 
the growth-rate is proportional to velocity gradients, the 
fastest growth is for the fluctuations shorter than the 
viscous scale of turbulence, they are the first to reach 
saturation and stroiigly infiuence the subsequent evolu- 
tion of the system ||7[]. It is then important to have a 
systematic description of the field that has emerged from 
the linear dynamo phase. Consistent description of the 
long-time evolution of the small-scale field with the ac- 
count of diffusion remained elusive for a long time [g|, 
only the second moment has been found g,g| . When dif- 
fusivity K is small, the field is almost frozen into the fiuid 
and is expected to grow exponentially like an infinitesi- 
mal material line element. To what extent this is offset 
by a transversal contraction that eventually brings diffu- 
sion into play depends on the statistics of stretching and 
contraction. We flnd below the growth rate of the field in 
a typical realization 7 = (log B)/t for an arbitrary veloc- 
ity statistics and derive analytically the whole function 
En — log (i3^")/t for a short-correlated velocity. Both £"„ 
and 7 are finite at k — > -|-0 (what is called fast dynamo 
p,B|). Particularly interesting is that 7 = Ai in a perfect 
conductor while 7 < Ai/2 at whatever small resistance, 
Ai being the largest Lyapunov exponent i.e. the growth 
rate of a material line element. 

Consider the kinematic stage of dynamo when the only 
equation to solve is that for the magnetic field 



diffusivity ratio is large and there is a wide interval of 
scales between viscous and diffusive cut-offs where the 
velocity is spatially smooth while the magnetic fleld has a 
nontrivial spatial structure to be described below. Given 
the initial condition, the solution of (nl) for the smooth 
case (when one may substitute v = or introducing the 
local strain aap = dva/drp) is most conveniently written 
in Fourier representation 



B(k, t) = W (t) B (k(0), 0) exp ~k / k'^{t') dt' 



dtB + (v • d)B = (B • a)v + kAB, 



(1) 



while the velocity statistics is presumed to be known. 
In many astrophysical applications the viscosity-to- 



where the wavevectors evolve as k(f') = W'^{t,t' )]<.{{) 
and the flnal condition is k(t) = k [nO[. The evolu- 
tion matrix W satisfles dW{t,t')/dt = a{t)W{t,t'), with 
W{t',t') = 1 and W{t) = W{t,0). We adopt here the 
methods of Lagrangian path analysis |l^,|ll] developed 
recently for the related problem of passive scalar ||l^,H . 
The moments of B are to be calculated by two inde- 
pendent averaging: flrst, (trivial) average over the ini- 
tial statistics and, second, average over velocity statis- 
tics. Without any loss of generality, we assume the ini- 
tial statistics to be homogeneous, isotropic and Gaussian, 
with zero mean and the variance (^^(k, 0)i?/3(k', 0)) = 
Pc./3(k)fc2 f{P)5(k + k'). The solenoidal projector P^f} = 
(pap — kakfj/k"^) ensures that B is divergence-free. The 
initial magnetic noise is concentrated at the scale L, we 
use /(A:^) = L^exp(— /c^L^) whenever an explicit calcu- 
lation is performed. Inhomogeneous advection produces 
smaller and smaller scales and balances with diffusion at 
the scale r^ = a/k/Ai, the magnetic Reynolds number 
L/rd is assumed to be large. 

The wavectors are of order 1/ L initially and for some 
period they remain much smaller than l/r^. This is 
the stage where the dynamics is insensitive to diffusion 
so that the fleld is frozen into the fluid like in a per- 
fect conductor, see [p|JTo[| and (||) below. At some time 
t(i oc ln(L/rrf), the wavevectors reach 1/r^, transversal 
contractions bring diffusion into play and the new regime 
starts which is the main subject of this paper. It is sup- 
posed that td, is much larger than the velocity correlation 



time r and we can then carry over from random matrix 
theory the well-known Iwasawa decomposition (see e.g. 
0) : W{t) = RDS. For any fixed time, R is an S0(3) 
rotation matrix, D is diagonal with T)ii(t) = exp[/9i(i)] 
and the shearing matrix S is upper triangular with unit 
elements on the diagonal. The sum of p^'s vanishes by 
incompressibility and the ratios pi/t tend at i ^ oo to 
the three Lyapunov exponents Xi (arranged in decreasing 
order) — see e.g. [Ts] ]. 

Since B is expressed via W, our aim now is to reformu- 
late the average over a into that over W. The average 
over R is equivalent to the integration over the direc- 
tions of the vectors involved. The matrix S tends for 
each realization to a time-independent form. Indeed, it 
will be shown below that the magnetic field moments 
grow exponentially in time and the dominant contribu- 
tions come from realizations with pi ^ /52 at i — > oo. For 
such realizations the matrix S is frozen at large times, the 
eigendirections of W^ W do not fluctuate in time and fix 
an orthogonal basis : rii , n2 , n^, with respective stretch- 
ing (contraction) rates 2Ai, 2A2, 2A3 [Q. We shall also 
see below that the time-independent random matrix ele- 
ments of S influence only constant factors, not of interest 
for the space-time dependencies studied here. The prob- 
lem is reduced then, on one hand, to the integration over 
the angles of the vectors involved and, on the other hand, 
to the average over the statistics of pi and p2 • 

The moments of the magnetic field can be obtained by 
considering B^(i) averaged over initial statistics 

then taking powers and averaging over the velocity. Here 
we have made the change of variables q = W^li that 
reexpresses the average in terms of the wavevectors at 
t — 0. When t 3> lAsj"^ the main contribution to 
A(i) = J^dt'W-\t')W-^^'^{t') is given by t' > \X3\-\ 

Changing the variables q = S''^Q we eliminate the con- 
stant S'-matrix from the diffusive exponent. The depen- 
dence on S remains just in the quadratic in Q pref- 
actor and in / so that in averaging over the velocity 
we may replace S by the unit matrix. It follows that 
qAq = /„* dt J2 Q?e-2p. = u{Q, p) and in the q^Tr term 
the main contribution is e^i^^{Q\ -t- Qg). In any given 
realization the growth of the field is thus described by a 
simple formula 

B2(i) ~ j d'QfiQ) exp{-2KU)e'''^ {QJ + QJ) . 

The dynamics of the field growth is now clearly seen. Ini- 
tially diffusion is unimportant {kU "C 1) and B^ grows 
as e'^''^ i.e. as a square of a material line element. At 
t ~ td = I A3 1 ~ log L/rrf, the diffusive exponent starts 
to decrease substantially and the growth rate is reduced. 
Asymptotically for i ^ t^, it is clear that the realiza- 
tions and the q's dominating the growth are such that the 



quadratic form kU at the diffusive exponent remain 0(1). 
Note that for growing functions one has with exponential 
accuracy/ (it'exp(— p3(t')) oc exp(— p3(i)). The integra- 
tion over Qs is thus restricted within exponentially small 
interval: the diffusion exponent remains 0(1) only for 
initial wavevectors with such a small projection on the 
contraction direction n3 that the respective component 
does not reach l/r^ during the time t |jl^. Neglecting 
Ql comparatively to Q2 and omitting numerical factors, 
we get after the integration over Q 

B2(i) ^ exp[2pi(t)]{l + (rd/L)2exph2p3(t)]}-'/' 

x\l + {rd/Lf I dt'eM~'^p2{t')]\ ■ (2) 



The first figure bracket reduces for large times to 
{L/rd) expp3, the geometrical factor due to the orthogo- 
nality condition to n3 . In the second line the exponential 
term can be either comparable or larger than unity de- 
pending on the sign of p2, which corresponds to the two 
geometrical pictures (cone vs pancake in /c-space) illus- 
trated in 0. 

Moments of (||) should be averaged over the prob- 
abihty distribution Vt{pi,P2)- When t » Aj"^,r, 
the theory of large deviations assures that Vt{pi, P2)o(- 
exp[— iiJ(pi/t, P2/O], where the entropy H{x,y) has a 
sharp minimum _ff = at a: = Ai, j/ = A2 whose width 
decreases as t~^^^ [Q (for the vector case see, e.g., Q). 
The mean growth rate ^{t) = (lni3^(t))/2t is then simply 
obtained by taking the logarithm of (H) and substituting 
there pi — Xit (strictly speaking, one had to average the 
logarithm over the initial measure as well, yet this dif- 
fers by a correction decreasing as i~^). At t <^ td the 
growth rate 7 = Ai, as it has to be for a perfect conduc- 
tor. During an intermediate stage t ^ td the growth rate 
decreases and eventually at i 3> id it comes to an asymp- 
totic value 7oo , independent of k (so called fast dynamo 



i{(Ai-A2)/2, (A2-A3)/2}, 



(3) 



Note that the growth-rate does not fluctuate at large 
times, 7 > and 700 — >• as A2 — > Ai or A2 -^ A3, 
corresponding to the zero growth rate for axially sym- 
metric cases. For a time-reversible flow statistics A2 = 
and 7oo = Ai/2. For isotropic Navier-Stokes turbulence, 
numerical data suggest A2 ~ Ai/4 |9|, so that our predic- 
tion for the long-time growth rate in a typical realization 
is 7oo — 3Ai/8. 

The moments with n > all grow in a random incom- 
pressible flow with a nonzero Lyapunov exponent since 
En ~ log(i3^")/2t is a convex function of n (due to 
Holder inequality) with Eq — and dEn / dn{0) = 7 > 0. 
Even when 7 = 0, the growth rates _E„ are positive for 
n > if -ff has a finite width, that is if the flow is ran- 
dom (for n = 1 this was previously stated in ||l^). The 



growth of the 2n-th moment at i ^ i^ is determined by 
the average of exp(2npi). For t ^ td, the expression to 
average is either exp(n/9i — np2) (with p2 > ^ogrd/L) 
or exp(np2 ~ i^Ps) depending on whether the entropy 
function favours positive or negative p2 (cf. |l^])- The 
function En is nonuniversal since it is determined by the 
saddle-point of 'Dpi,2 integration which depends on the 
particular form of the entropy function. The saddle-point 
falls within the (universal) parabolic region of H around 
the minimum only for n <C (Ait)~^. Therefore, the tem- 
poral growth of the moments will be calculated below for 
a short-correlated velocity. 

Here, we continue with the general case to establish 
what is universal in the different-point correlation func- 
tions F2n = (niLi (B(x2fe„i,t)B(x2fe,t))). Its calcula- 
tion is reduced to averaging (2n — 1)!! terms arising from 
the Wick decomposition in the Gaussian integration over 
the random initial condition. Each term is a product of 
n integrals generalizing the one for B^ with the inclusion 
of respective exp(zrjV[^-^'~^q) in the integrand. The n 
vectors r^ are differences between couples of x^'s. The 
new feature with respect to the moments is the presence 
of the rotation matrix R in the exponential factor. The 
q-integrations proceed along the same lines as previously: 
in any of the n integrals we change variables q = iS-^Q 
and the dependence on S is thus entirely moved into 
the prefactors. Substituting S by the unit matrix and 
performing the Q-integrations, we obtain the long-time 
asymptotics for any of the (2n— 1)!! contributions to F2n '■ 



(L/rd)" exp [7i(pi - P2)] 



-2ps 



n 



2- 



rj/L'- 
L2 



i5n/2 



dcosi 



exp 




Rlj 



df / d(j) (4) 
"'0 



Rlj 



4 [L2e2p2 + 2rj] Sr^ 



where R^- — R3[(f]R2[6]R3[(j)]rj, and ^2,3 stand for ro- 
tation around Y and Z axis respectively. We consider 
time-reversible statistics, where the scaling laws turn out 
to be universal. Let us explain the physical meaning 
of (0) and derive the correlation functions starting from 
n — I. The realizations contributing have the advec- 
tive exponent exTp{irW^'~^q) of order unity. This re- 
quires p2 > ln(r/L) and the direction of contraction n^ 
to be almost perpendicular to r, which gives angular fac- 
tor (L/r)expp3. At Xit > \n{L/r) we get then 



F2(r,t)~- / dpi / dp: 

^ J-00 J\n(r/L) 



.^ePi~P2~tH ^^-2-h^E,t^ 



where h= dH/dy taken at y = and at x given by 
the saddle point dH/dx = 1. Time reversibility means 
that H{x,y) = H{x + y,—y) so h = 1/2 and F2 ex 
r~^/^eyi-p{E2t). At r < L and Aif-C ln(L/r), F2 is r- 
independent. This generalizes the consideration of M for 
an arbitrary time-reversible statistics. To understand the 



simple geometrical picture behind this derivation, note 
that the integral over p2 comes from p2 — ln(r/L) . That 
means that the field configurations in the form of strips 
with width r dominate F2{r). Angular integral in (Q) 
comes from (^ ~ 1 that is the strips with the stretching 
direction almost parallel to r do not contribute (because 
of cancellations due to solenoidality) . 

For n > 2, the geometry of the vectors Vj becomes im- 
portant. Let us first consider the case where all the vec- 
tors are in the same plane (their length being r). They 
can be either coUinear or not. Almost orthogonality of 
Tj to na involves therefore either one angle or two, giving 
the angular factor {L/r) exp p^ or its square, respectively. 
The other difference concerns the behavior along n2 . For 
non-coUinear geometry not all vectors can be orthogonal 
to n2 and p2 should then be constrained as p2 > log r/L. 
This is technically signalled by the fact that the integra- 
tion over if is not saddle-point. Conversely, for coUinear 
geometry all the vectors can be orthogonal both to n2 
and n3 , giving an additional angular factor (L/r) exp p2 : 
the saddle-point integrations over 6 and ip pick 6 — 7r/2, 
ip ~ eKp[p2]L/r. In the rest of the integrals (either n ~ 1 
or n — 2, respectively) the above angular constraints en- 
sure that the advective exponents are 0(1) so the dif- 
fusive exponents exp(— 2KqAq) become important. The 
calculation of these integrals is essentially the same as for 
the moments and this is where diffusion comes into play. 
The wavevectors should be quasi orthogonal to n^ , giving 
either n — 1 or n — 2 factors (L/rd) exp p^. The growth 
along n2 for a generic planar geometry is automatically 
controlled by the previous constraint p2 > log r/L; for 
coUinear geometry it provides the bound p2 > logrd/L. 
Simply speaking, the strips with the width rj stretched 
along r contribute in the coUinear case, while the width 
is r in the generic case. The resulting integrations over 
pi are saddle-point and those over p2 are dominated by 
the lower bounds. Finally 



F2n 



„E^t 



1 

(rd/r)' 



coUinear 
2 planar. 



(5) 



Here, £'„ = a;„ — H{xn,0) with dH/dx{xn,Q) = n. The 
fact that the integrals over p2 are all dominated by the 
lower bounds indicates the geometric nature of the scal- 
ing universality found : the magnetic field configurations 
that contribute are narrow strips (not ropes and layers 
suggested in [nO[) with one direction of stretching, one of 
contraction and a neutral one. The factor {rd/r)"^ is the 
probability for two points at distance r to lie within the 
same strip of width rd ■ 

The peculiar nature of strip configurations has another 
dramatic consequence for n > 3: the correlation func- 
tions are strongly suppressed in a generic situation when 
at least three vectors r^'s do not lie in parallel planes. In- 
deed, they cannot then be on parallel strips and nonzero 
correlation appears only because the strips have expo- 



nential diffusive tails. As a result, the factor (rd/r)"^ in 
the planar formula is replaced by exp(— ar^sin 0/r^), 
where a ~ 1 and is the minimal angle between a vec- 
tor and the plane formed by another two. That can 
be derived from (0) where all angular integrations are 
not saddle-point now. For a general irreversible velocity 
statistics, the r-dependencies are different yet the quali- 
tative conclusions (that the correlation functions are not 
exponentially suppressed only for planar geometry and 
are anomalously large for coUinear geometry) are gen- 
erally valid. Note that angular anomalies are likely to 
be peculiar of the viscous interval where advection by a 
smooth velocity preserves coUinearity. Similar coUinear 
anomalies (which are manifestations of the strip con- 
figurations described here) have been found before for 
a passive scalar advected by a smooth velocity pO| , pi| . 
Note that the cliff-and-ramp structures observed in pas- 
sive scalar experiments (see [^ for review) are probably 
related to the strips too. 

It is left to find En for the standard Kazantsev- 
Kraichnan model of an isotropic short-correlated Gaus- 
sian strain with (crQ/3(i)cra/3(0)) ~ 10Xi6(t). Straight- 
forward derivation gives Gaussian T't{pi,P2) with H ~ 
(pi -f- P2/2 - Xitf/Xit + ipl/AXit. Note that A2 = 0. 
Now we integrate the moments of (g) with such Vt- First, 
we integrate exp(2npi) and get the answer for the perfect 
conductor pjl0|| 



(B^")- exp [2Xin{2n + 3)i/3] 



(6) 



The main contribution to the n-th moment comes from 
pi — Aii(4n -|- 3)/3,p2 — — 2nAii/3 so (||) is vahd un- 
til Lexppa > Td that is for t < Md/{2n -f- 3) with 
td = A^"'^ \n{L/rd). Then there is a logarithmically wide 
cross-over interval when the growth is non-exponential. 
The asymptotic regime starts at i > 3td/{n + 2) when 
unity in the first parenthesis of (|^) may be neglected. The 
integral over pi now comes from pi +P2/2 = (n + 2)Xit/2 
while that over p2 is dominated by p2 ~ In^rd/L): 

(B2") ~ (i/r,)^"/'exp(£;„t), E„ = Ain(n + 4)i/4 . (7) 

For n — 1, this was obtained by Kazantsev ||]. The 
difference between (^) and (0) formally means that the 
two limits t — > cxD and k ^ do not commute (what 
is called dissipative anomaly). The physical reason is 
quite clear : realizations with continuing contraction 
along two directions contribute most in a perfect con- 
ductor, while with diffusion present one direction is neu- 
tral. A magnetic field initially concentrated in the ball 
with the radius L will have the fastest growth rate 
7 if the ball turns into a strip with the dimensions 
i(L/rd)i/2exp(3Di), rd, L{rd/ Lf'^ c^-iDt). The 



main contribution into the n-th moment will be given by 
strips with the dimensions L[LlrdS^I'^ exp[3(n-|-2)Di/2], 
Td and L{rd/LY''^ exp[-3(n + 2)Di/2]. To conclude, we 
related the growth-rate of the small-scale dynamo to the 
Lypunov exponents of the flow and described analytically 
the strip structure of magnetic field. 
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